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PLANAR HETEROSTRUCTURE 
GRAPHENE — NARROW-GAP SEMICONDUCTOR — GRAPHENE 

P. V. Ratnikov^ and A. P. Silin 

We investigate a planar hetero structure composed of two graphene films 
separated by a narrow-gap semiconductor ribbon. We show that there 
is no the Klein paradox when the Dirac points of the Brillouin zone 
of graphene are in a band gap of a narrow-gap semiconductor. There 
is the energy range depending on an angle of incidence, in which the 
above-barrier damped solution exists. Therefore, this heterostructure is 
a 'filter" transmitting particles in a certain range of angles of incidence 
upon a potential barrier. We discuss the possibility of application of this 
heterostructure as a "switch". 

Graphene is a two-dimensional gapless semiconductor, and charge carriers are massless 
Dirac fermions [1]. It is known [2] that a massless relativistic particle with spin 1/2 possesses 
the chirality property, i.e. it is characterized by a certain spin projection onto its momentum. 
In case of graphene, the chirality is defined by a projection of pseudospin onto a momentum 
direction, which is positive for electrons but negative for holes near K point of the Brillouin 
zone (BZ) [3], i.e. an electron and a hole are analogues of a massless neutrino with the right- 
and left-hand helicities, respectively. However, the situation is inverse near K' point where 
electrons and holes have the left- and right-hand helicities, respectively, [4, 5]. The massless 
relativistic particle is described by one spinor, i.e. a two-component wave function [6, 7]. It 
gives basis to state that the effective Hamiltonian describing the charge carriers in graphene 
near K point is a 2 x 2 matrix, and corresponding equation is the Weyl equation^ 

ucr ■ pip = Eip, (1) 

where u = 9.84 x 10^ cm/s is the Fermi velocity which is an analogue of the Kane matrix 
element for the rate of interband transitions in the Dirac model [10], p = — i^V (hereafter 

^ r at nikov@ lpi.ru 

^The massless fermions considered separately near K and K' points are similar to the Weyl (two-component) 
neutrino. The Dirac equation is used for its simultaneous description [8]. The Dirac equation is equivalent to 
a pair of the Weyl equations. The Dirac equation in two-dimensions can be written as a 2 x 2 matrix (it can 
be used equivalently with a 4 x 4 matrix representation) which coincides with the Weyl equation for a massless 
particle on a plane. However, the former equation can be also used for a description of a particle with finite 
mass. Using this fact, the problem very close to considered in this paper task was earlier solved by Gomes and 
Peres [9]. 
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^ = 1), and cr = (cr^, ay) are the Pauli matrices. The dispersion relation of the charge carriers 
is hnear in momentum k 

E = ±uk. (2) 
A narrow-gap semiconductor is described by the 4x4 matrix Dirac equation [11] 

Hd"^ = {mck ■ p + /3A + Vo} ^ = E^, (3) 
where \1/ is a bispinor, u is the Kane matrix element for the rate of interband transitions, 




^0 -/^ 

the band gap, and Vq is the difference of work functions of the narrow- gap semiconductor and 
graphene (|Vo| < A). 

It is necessary to introduce the four-component wave function, bispinor, for simultaneous 
description of the charge carriers in graphene and the narrow-gap semiconductor. In this case, 
the Dirac Hamiltonian is 

So^( V "^P). (4) 

Hamiltonian (4) is equivalent to Hamiltonian used in Ref. [8], with an accuracy of two consec- 
utively performed unitary transformations U2 = -i) Ui = (0^!,) [12]. In addition, 
both left-hand and right-hand helicity massless fermions are presented in the system. Tran- 
sitions between K and K' points are improbable, so it is possible to consider that particles 
conserve the chirality property. 

By performing the unitary transformation U2-, it is convenient to present the Dirac Hamilto- 
nian describing the charge carriers in the total heterostructure (Fig. 1(a)) in the form in which 
diagonal blocks contain momentum operators 

H'^=h--'^ + ^^ ^1 V (5) 

y Ai -UiO- -"p + Vij 

where Ui = = u, Vi = V3 = 0, and Ai = A3 = are the parameters related to graphene, 
U2 = u, V2 = Vq, and A2 = A are the parameters of the narrow-gap semiconductor (Fig. 1(b)). 

For the components of the bispinor describing a particle in graphene, the following equalities 
exist 

^2 = s^ie^^ 

(61 

^4 = -#3e'*, 

where = arctan is the polar angle of momentum k = (kx, ky) of the charge carriers in 
graphene (the angle of incidence), s = signE. 

For the components of the bispinor describing a particle in the narrow-gap semiconductor, 
the following equalities exist 

, E - Vq ^ uqx- iuky , 

V'S = T T ^2, 

_^ _ ^ (7) 

UQx + luky E -Vq 
i'i = ^ ^1 + ^2, 
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where 

u'ql = {E- V,f -A'- uX- (8) 

We find the solution in three ranges: I) x < 0, II) < x < III) x > D {D the width 
of the narrow-gap semiconductor ribbon, see Fig. 1(a)), taking into account relations (6), (7) 
and assuming that the solution is oscillating in range II (g^ > 0), 
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Fig. 1. Considered planar heterostructure; (a) two graphene layers separated by the narrow-gap 
semiconductor ribbon with the width D (it is hatched), (b) the band structure: level E = 
corresponds to the position of the Dirac points in the BZ of graphene, the band gap of the 
narrow-gap semiconductor is Eg = 2A, Vq is the difference of the work functions of graphene 
and the narrow-gap semiconductor, completely filled valence bands are hatched. 
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where r and t are the reflection and transmission coefficients, respectively, [3], ci, C2, oi, 
02, bi, 62 are complex constants determined from the boundary conditions'^. 
Using the boundary conditions [13, 14] 

v^^(-) = V^¥+\ (12) 

where the quantities marked by "(— )" and "(+)" relate to material placing on the left and 
right of the boundary, respectively, we obtain for the transmission coefficient 



COS cos{qxD) + i ^tan 6* sin — s ^-^^ j sin^q^D) 



^-^k.U^ (13) 



where tan^ = Expression (13) corresponds to the oscillating solution in range II. In 
order to get the transmission coefficient at the exponentially damped solution in range II, the 
replacement qx — ^ iqx should be made where v?q^ = + v?ky — {E — Vq)^, and q^ > 0. The 
transmission probabilities T = |tp for both kinds of solution in range II are 



cos^ 



Hascil — \ 2 ' i^^) 

cos^ 0cos^(g2:D) + ( tan 9 sin (p — s ^S^'^ j sm'^{qxD) 



cos^ 



^damp — , N 2 • i^^) 

cos20 cosh^ (qxD) + (^siiKp - s^;^] sinh^iqxD) 

One can see from formula (14) that Toscii = 1 when qxD = ttN, where is integer. It 
corresponds to maxima of the transmission probability shown in Fig. 2 (a)-(d). 

As one should expect, the transmission probability in case of the damped solution in range 
II is exponentially small for sufficiently large width of the narrow-gap semiconductor ribbon 
D > l/\qx\: Tdamp ~ e-^l?-!^. The result of the passage to limit, A ^ 0, in (13) coincides with 
the transmission coefficient t in Ref. [3]. 

The reflection coefficient is simply obtained 



-i sin{qxD) 



cos(0 -9)- e-'^ + ge^'^ ^--^°+^ 



cos cos 9 cos{qxD) + i (sin sin 61 - s^^) sin(g^L>) e"*^ - se-^<^^^=^|^ ' 

(16) 

where k' = a/^1~+~^- The passage to limit, A —>■ 0, in (16) is performed by replacements 
^ ^~^k^^ ^ = sign{E - Vq), the result coincides with formula (4) of Ref. [3]. 
The reflection probabilities R = |rp for both types of solution in range II are 



•^It should be emphasized that C2 = in the bispinors and for the right-hand hehcity particle, since, 
the equality ^^—p^'^a — "^r, where 75 — i/3, is valid for the bispinor describing the right-hand helicity 
particle, and Ci — and the equality ^+^ts vj/^ = is valid in those bispinors for the left-hand helicity particle 
[2]. Consequently, corresponding components of ^'7/ are zero on interfaces (ones are zero everywhere for the 
damped solution). 
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r, [cOs(0-g)-S^]^ 

cos^ cos^ 9 coi {qxD) + (sin sin 6 — s— ^) 
1 + 2s^^ cosU + 9)+ (^^^ 



y, MKT \-r ' I • M-'fc'^ /-| y\ 



_ n^gg cos^ + {uky sin - s(E - Vq))^ 

M cos^ cotli [fixU] + (w/Cy sin0 — — Voj) 
One can see from (18) that -Rdamp — * 1 at Ig^l-D ^ 1. It is simply verified that the following 
equalities are valid 



-^oscil ~l~ -Roscil 1 



(19) 

-^damp ~l~ -'^damp !• 




Fig. 2. The dependence of the probability Tosdi of the electron transmission through rectangular 
barrier being the band gap of the narrow-gap semiconductor GaAs with A = 705 meV on the 
angle of incidence, u = \J~-^ = 1-35 x 10® cm/s where m* = 0.068mo,mo is the free electron 
mass [16]. The difference of work functions of GaAs and graphene is assumed to be positive 
and equal in Vq = 100 meV. The angle 0o ~ 46.8° corresponding to the equality sin0o = u/u is 
marked. Two values of energy satisfying the above-barrier transmission condition > A + Vq 
are considered. When the angle of incidence approaches 0i, the upper boundary of the above- 
harrier damped range comes up to the energy E of a incident electron, for E = 1 eF 0i ~ 24°, 
for E = 2 eV <pi ^ 40°.- (a) E = 1 eV, D = 50 A; (b) E = 1 eV, D = 60 A; (c) E = 2 eV, 
D = 50 A; (d) E = 2 eV, D = 60 A. 
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Let us analyze the conditions at which the oscillating or damped solution can exist in range 
II. For definiteness, we consider the case of electrons^: its energy E = uk is positive in graphene, 
then, for the oscillating solution, the following equality should hold: 



uk = Vo + + + u^k^ (20) 

which is valid at condition 

uk-Vo> ^J^^+lKl- (21) 
Conversely, it is necessary for the damped solution^ that 



uk = VQ + y^A2 - vi^ql + m2^2_ (22) 

It is valid at condition of the intersection of the dispersion curves for graphene and the narrow- 
gap semiconductor [15] 

uk -Vq< ^^^^u^kl. (23) 

It is evident from the inequality (23) that if the Dirac point of the BZ of graphene falls into 
the band gap of the narrow-gap semiconductor (tunneling through the potential barrier being 
the band gap of the narrow-gap semiconductor) then the solution in range II for electrons with 
the energy i?e < Vq + A (analogously for holes with the energy i^/i > Vq — A) is always damped 
one. 

The momentum range corresponding to the oscillating solution is defined by inequality 

{u^ - sin^ 0) k^ - 2Vouk + - A^ > 0, (24) 

and the inverse inequality defines the momentum range of the damped solution. The analysis 
of inequality (24) shows the following: 

1) if M > M then at any angle of incidence — | < < | for electrons with energy Eg and 
holes with energy Eh in range 

A + Vo<E,<E+{c^), 

(25) 

E-{(P) <Eh< -A + Vo, 

where E^[(j)) = i-'!; s'iri^^^^ ^' ' ~ t?' ^^ere is the above-barrier damped solution; in 

range E^ > E^{(j)) and Eh < E^{(j)) there exists the oscillating solution; 

2) if u <u (it is valid for a number of the narrow-gap semiconductors, e.g., GaAs and InSb) 
then it is necessary to distinguish the following particular cases: 

a) the situation in the range of angles | sin0| < = is the same as the case 1); 

^Thc case of holes is equivalent to the case of electrons with an accuracy of the replacement E — > —E and 

Vo ^ -Vo. 

■^The expression with sign minus before the square root can also correspond to the damped solution if Vq > 
and value of the square root is smaller Vq. 
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b) the behavior of particles in the range of angles = < | sin0| < 1 is various depending 
on the value Vq: 



b*) if Ayl — |j < |Vo| < A then it should distinguish the subcases for all values of 
angles from this range: 

(i) there is the damped solution at any k (at any energy) for electrons at Vq > and 
for holes at Vq < 0; 

(ii) there is energy range above the barrier for electrons at Vq < and holes at 
Vo > 0, transparency window, in which there is the oscillating solution, and beyond it there is 
the damped solution: 

E,{<P)<E,<E,{<P), ^^^^ 
-^2(0) <Eh< -Ei{<P), 

where i^i ^(0) = 



u 


A 


"V 
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b**) if < A^l - p then 

(j) the situation is the same as case b*) in the range of angles = < | sin 01 < 

(jj) there exists only the damped solution at any k in range of angles " 

V " '0 

I sin 01 < 1. 

The potential barrier is an ideal reflector at sufficiently large angles of incidence in cases (i) 
and (jj), i.e. an "angle filter" transmitting particles with angles of incidence near to = 0. At 
the same time, it is supposed that \qx\E) ^ 1, i.e. T^amp ^ 1- Such an unusual feature of the 
rectangular potential barrier is related to the circumstance that the "speed of light" , analogues 
of which are u and m, is different in graphene and the narrow-gap semiconductor [14]. 

The case u = u should be attributed to the case 1). Then the energy range of the above- 
barrier damped solution disappears, and a particle behaves as an usual nonrelativistic particle, 
namely, there are the damped and oscillating solutions under and above the barrier, respectively. 
Similar results in this particular case have been obtained by Gomes and Peres [9]. 

Let us consider separately the case when we have instead of the narrow-gap semiconductor, 
a gapless semiconductor for which u u, and Vq 7^ (at m = m this case coincides with one 
considered in [3]). However, in contrast to Ref. [3], there is a number of features distinguished 
the case u ^ u. The transmission probabilities for both types of solution in the gapless semi- 
conductor are given by expressions (14) and (15), the only difference is that it is necessary 
to make the replacement E — Vq ^ suk' . In the above manner, let us analyze what kind of 
solution we have in the gapless semiconductor: 

1) if M > M then at any angle — | < < | 

a) there exists the oscillating solution for electrons at Vq < and for holes at Vq > for 

any fc; 

b) there exists the damped solution for electrons at Vq > and for holes at Vq < in 
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the energy intervals 

E^{cl>)<E,<E,{ct>), ^^^^ 
E, {({>)< Eh <E+ 

where E^{(j)) = Vq; and there exists the oscillating solution beyond these intervals. 

If we regard Vq as the potential barrier height [3], then we have the under-barrier oscillating 
solution, this fact corresponds to the Klein paradox; 
2) if u <u then 

a) the situation is the same as in 1) for angles | sin0| < |; 

b) for angles | < | sin0| < 1, we are should distinguish two particular cases: 

(i) the solution is damped for electrons at Vq < and holes at Vq > for any k; 

(ii) the solution is oscillating for electrons at Vq > and holes at Vq < in energy 
ranges (27) but out of ones there is the damped solution. 

Finally, let us consider the particular case A = and Vq = at u ^ u: 

1) if u > u then the solution is oscillating at any angle — | < < | and any energy, this 
fact corresponds to the Klein paradox; 

2) if u < u then the solution is oscillating at any energy for | sin0| < | and the solution is 
damped at | sin0| > |. 

In conclusion, we note that the considered heterostructure can be used as the "switch", 
namely, applying a voltage on the narrow-gap semiconductor ribbon we can "switch on" and 
"switch off" transmission of the charge carriers through range II depending on the energy 
range in which the Dirac point of graphene falls (in the range of the oscillating or damped 
solution). When we apply an electric field F, the Dirac point of graphene shifts in energy 
by the value ~ eFd where is a distance from the voltage applying point to the narrow-gap 
semiconductor ribbon. We suppose that the electric field is weak enough: eFd < A — |Vb|, 
i.e. current does not flow at the given Vq. The electric field correction results in displacement 
~ \eFD of extrema of the conduction and valence bands of the narrow-gap semiconductor 
[17]. Applying the voltage —Uq to the narrow-gap semiconductor ribbon changes the difference 
Vq = Vq — Uq of work functions between the narrow-gap semiconductor and graphene so that 
passage of electrons becomes possible at eFd > E^{(l))\v;^. Condition of passage for holes is 
eFd > \E^{(j))\. Changing Uq, we can achieve passage of either electrons or holes. 

An alternative scheme of the "switch" is possible. Due to zero gap in graphene one can pump 
electrons from the substrate in the conduction band or displace electrons from graphene thereby 
obtaining holes in the valence band. Changing position of the Fermi level Ep in one of the 
graphene layers we can provide passage of either electrons at the condition eFd + Ep > E'^{(f)) 
{Ep > 0) or holes at -eFd + Ep < E-{(f)) {Ep < 0). 



8 



References 

[1] K. S. Novoselov et al., Nature 438, 197 (2005). 

[2] S. S. Schweber, Introduction to Relativistic Quantum Field Theory (Harper and Row, 
New York, 1961). 

[3] M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nature Physics 2, 620 (2006). 
[4] P. L. McEuen et al., Phys. Rev. Lett. 83, 5098 (1999). 

[5] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, 
arXiv:0709.1163v2. 

[6] L. D. Landau, ZliETF 32, 405 (1957) [Sov. Phys. JETP 5, 336 (1957)]; 
L. D. Landau, Nucl. Phys. 3, 127 (1957). 

[7] A. Salam, Nuovo Cim. 5, 299 (1957). 

[8] T. Ando, J. Phys. Soc. Japan 74, 777 (2005). 

[9] J. V. Gomes and N. M. R. Peres, J. Phys.: Condens. Matter 20, 325221 (2008). 

[10] P. V. Ratnikov and A. P. Silin, Kratk. Soobshch. Fiz. No. 1, 46 (2008) [Bull. Lebedev 
Phys. Inst. 35, 27 (2008)]; arXiv:0801.4642vl. 

[11] B. A. Volkov, B. G. Idlis, and M. Sh. Usmanov, Usp. Fiz. Nauk 65, 799 (1995) [Phys. Usp. 
38, 761 (1995)]. 

[12] P. V. Ratnikov, Pisma v ZhETF 87, 343 (2008) [JETP Lett. 87, 292 (2008)]; 
arXiv:0805.4510vl. 

[13] A. P. Silin and S. V. Shubenkov, Fiz. Tverd. Tela 40, 1345 (1998) [Phys. Sol. St. 40, 1223 
(1998)]. 

[14] A. V. Kolesnikov and A. P. Silin, ZhETF 109, 2125 (1996) [JETP 82, 1145 (1996)]. 

[15] A. V. Kolesnikov, R. Lipperheide, A. P. Silin, and V. Wille, Europhys. Letters 43, 331 
(1998). 

[16] A. P. Silin and S. V. Shubenkov, Kratk. Soobshch. Fiz. No. 7-8, 9 (1996) [Bull. Lebedev 
Phys. Inst. 23, 8 (1996)]. 

[17] P. V. Ratnikov and A. P. Silin, Kratk. Soobshch. Fiz. No. 11, 22 (2005) [Bull. Lebedev 
Phys. Inst. 32, 17 (2005)]. 

Received 29 May 2008. 



9 



